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In the high temperature, chirally invariant phase of QCD, 
the quark propagator is shown to have two sets of poles with 
different dispersion relations. A reflection property in mo- 
mentum space relates all derivatives at zero-momentum of 
the particle and hole energies, the particle and hole damping 
rates, and the particle and hole residues. No use is made of 
perturbation theory. 
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fermions [gSj. Farrar and Shaposhnikov exploited the 
importance of chirahty transport in baryogenesis at the 
electroweak scale and used the doubled fermion disper- 
sion relations to account for the baryon asymmetry of the 
universe |2^. However, it has been shown that includ- 
ing fermion damping makes the effect too small [^0|-^. 
There are renewed attempts to account for baryogene- 
sis at the electroweak scale using the fermion dispersion 
relations in models with two Higgs doublets |34| and in 
minimal supersymmetric models Bq). 



1. Background 



2. QCD Motivation 
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In the high temperature, chirally invariant phase of 
gauge theories the fermion propagator has some unusual 
properties in the one-loop approximation. Despite ex- 
plicit chiral invariance the fermion has an effective mass 
proportional to temperature [^^J . The "mass" is a con- 
sequence of the preferred reference frame of the heat 
bath, which allows the pole of the fermion propagator 
to be off the light cone. The second and more surpris- 
ing effect is that there are two poles in the propagator 
corresponding to two different dispersion relations, both 
with positive energy . At zero momentum the two 

branches coincide. The branch describing normal parti- 
cle excitations increases monotonically with momentum. 
The other branch is a collective excitation, referred to 
either as a hole or as a plasmino, that decreases slightly 
at small momentum, reaches an absolute minimum, and 
then rises. The excitations on the upper branch of the 
dispersion curve have the same chirality and helicity (i.e. 
both positive or both negative) as is customary for par- 
ticles. The hole excitations of the lower branch have chi- 
rality opposite to helicity. These results are reviewed by 
Le Bellac §]. 

In the hard thermal loop approximation of Braaten and 
Pisarski, the effective fermion propagator is the starting 
point for consistent higher order calculations . The 

doubled fermion dispersion relations have been impor- 
tant in calculations of dilcpton production ]l^-p^ and 
strangeness production |15|| by a quark-gluon plasma. 
There have been investigations of how the dispersion re- 
lations are affected by retaining non-leading powers of 
temperature p6| , p^ , by including a bare mass p8| , pj[ , 
and by a including a chemical potential |20-22|. 

In the electroweak sector at high temperature the 
quarks, charged leptons, and neutrinos all have effec- 
tive thermal masses and doubled dispersion relations 
|^,^-27|. The phenomena also occurs with Majorana 



The evidence discussed thus far for the thermal fermion 
mass and for the doubled dispersion relations is based en- 
tirely on one-loop perturbation theory. The one-loop cal- 
culations are sufficiently accurate for electroweak effects 
but not for QCD phenomenology. Existing calculations 
are only valid if g is very small. Unfortunately, the most 
interesting aspects of the quark dispersion relations are 
limited both by the small coupling and also by some ac- 
cidentally small coefficients. For example, the thermal 
mass is TO = QA\gT . The minimum in the dispersion 
relation for the hole excitation is slightly below this at 
E = O.SSgT and the minimum occurs at a very small 
momentum, p = Q.\7gT. At temperatures well above 
the critical temperature g is indeed small, but near it 
is not. 

There is some nonperturbative evidence that the one- 
loop calculations are qualitatively correct. Peshier et al 
p6| , ^ and also Levai and Heinz |Q have successfully 
fit lattice simulations of high temperature QCD using 
effective quark masses [|3^-|3^. These studies do not 
test the doubled dispersion relation. Recently Schafer 
and Thoma have computed the qua rk self-energy in the 
presence of a gluon condensate p}9[| . Their calculation 
is one-loop but with the gluon condensate determining 
the gluon propagator. In their calculation it is the soft 
loop momenta that control the quark dispersion relation 
rather than the hard momenta as is usually the case. 
Nevertheless the results are qualitatively similar in sev- 
eral respects. There is an effective mass, but now related 
to the condensates {E"^) and (S^). Using lattice data for 
the condensates gives m « LIST for temperatures in the 
range LlTc < T < ATc- There are, in addition, two dif- 
ferent dispersion relations. They coincide at p = but at 
small p have opposite slopes. The dispersion relation for 
the quark is monotonically increasing; that for the hole 
decreases slightly, reaching a minimum at p ~ to/ 2 and 
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then increasing at larger p. The self-energy computed by 
Schafer and Thoma is a polynomial in pq and p divided 
by {pq — p^Y, whereas the hard thermal loop self-energy 
contains ln(|)o ip)- 

This naturally leads to the question of whether the 
existence of the separate particle and hole dispersion re- 
lation are independent of weak coupling and of the one- 
loop approximation. This paper will demonstrate that, 
independently of perturbation theory, the fermion prop- 
agator will always have two distinct dispersion relations 
in the high temperature, chirally symmetric phase. 



3. Retarded propagator 

It will be most convenient to analyze the retarded prop- 
agator and later express the time-ordered propagator in 
terms of it. The free retarded propagator is 



loPo -1 -P 



{Po + iriY 



and is independent of temperature. It can be rewritten 
as 



ShcciPO,P) 



l7o-7-p l7o+7-P 



2po-p + ij] 2po+p + iri 



(1) 



The particle pole is in the fourth quadrant of the complex 
Po plane at po — p—if]; the antiparticle pole is in the third 
quadrant at pa — — p — irj. 

The full retarded propagator has both poles and cuts 
but the singularities are always in the lower half-plane. 
At nonzero temperature the propagator in the rest frame 
of the heat bath will depend separately on energy po and 
momentum p. Invariance under chirality and parity lim- 
its the self-energy to be a linear combination of the ma- 
trices 7o and 7 • p. The most general possibility can be 
written 



D+{po,p) 



1 -P) ^ ^(70 + 7 -P) 
D-{pa,p) 



(2) 



The arguments below will show that contains four 
poles. In the fourth quadrant there is a pole at a com- 
plex energy £p for particle excitations and another at a 
different complex energy £h for hole excitations. The two 
poles for the corresponding antiparticle and antihole are 
in the third quadrant. 



4- Existence of separate dispersion relations 

The first step in the argument is the observation that 
at p = the propagator cannot depend on the direction 
of the unit vector p. This is because the dependence on 
p of the fermion self-energy comes entirely from the the 
vector p as shown in Appendix A. Therefore 13+ (po, 0) — 
D-{po,Q) and the propagator has the form 



S'^iPo^O)^ 



7o 



7o 



2D+{po,0) 2i?_(po,0) 



(3) 



Now comes the one dynamical input, namely that when 
p = there is a pole in the energy variable at the thermal 
mass of the quark so that 
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2{po-M)a{po)' 



(4) 



where a(po) is analytic at pq — M. The mass M is com- 
plex with a negative imaginary part since it is a pole of 
the retarded propagator. Denote this mass by 



A4 = m — 17/2. 



(5) 



From dimensional analysis both m and 7 must be pro- 
portional to temperature. The real part, m, has been 
extracted from lattice calculations pl[ |. 

The next step is to examine p^O. The denominators 
D±{pa,p) must have the structure 



D+{po,p) 
D-{po,p) 



{Po 



-M)a{pQ) 
-M)a{po) 



b{Po,p) 
c(po,p), 



where b and c both vanish at p = 0. Since b and c are 
continuous in p and vanish at zero momentum, they can 
be made arbitrarily small by choosing p sufficiently small. 
Therefore in the neighborhood of po ~ Al and p « one 
can approximate 



D+ipQ,p) 
D-{po,p) 



{Po 
{Po 



M)a{M) 
M)aM) - 



'b{M,p) 
c{M,p). 



vanishes at a complex energy po = £p given by 
b{M,p) 



£p = M 



a{M) 



+ . . . (p small) . 



(6a) 



The subscript on £p indicates that it is the complex en- 
ergy of the particle excitation. It is, of course, a function 
of the momentum p. D_ vanishes at a different energy 
Po = £h given by 



Sh = M 



c{M,p) 
aiM) 



(p small). 



(6b) 



which is the complex energy of the collective hole exci- 
tation. Note that it is not necessary to assume that b 
and c are differentiable at p = 0. Equations (^^ and 
( pq ) must have different momentum dependence because 
b(M , p) 7^ c{M , p) or else the full propagator in Eq. (||) 
would be completely independent of 7 • p. 

For general momentum p the particle and hole energies 
are solutions to 



M)ai£p) 
-M)ai£h) 



b{£p,p) 
c{£h,p)- 



(7a) 
(7b) 



Since these energies are poles of the retarded propagator, 
they have negative imaginary parts: 
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Op 



Epip)-ijp{p)/2 
Eh{p)~iih{p)/2. 



(8a) 
(8b) 



At zero momentum i?p(0) = i?/i(0) and 7p(0) = 7/1 (0) as 
indicated in Eq. (^). Although the arguments for this 
depend only on rotational invariance, perturbative cal- 
culations do show these properties. Equality of the zero- 
momentum energies is displayed in the one-loop calcula- 
tions [p]J^-|^,^ ; equality of the zero- momentum damping 
rates was found by Braaten and Pisarski ||^ using HTL 
resummation. 



5. Reflection symmetry 

The momentum variable p originally has a positive, 
real value \p\. However, one can extend the functions 
D±iPOiP) to more general values oi p. In particular, it 
is useful to allow p to be real and negative. Appendix B 
shows that 



D-{po,p) = D+{po,-p), 



(9) 



or equivalently 6(poi ~p) = c{po,p). Under the change p 
to —p in Eq. (|7^), let the solution be = Ep{—p). Then 
Eq. (|7a|) becomes 



0= {VL- M)a{n)^c{yL,p). 



Comparison with Eq. (7b) shows that f2 
gives the reflection relation 

The real and imaginary parts of this are 



£h{p). This 



(10) 



Ep{-P) 
Ip{-P) 



Eh[p) 
lh{p)- 



(11a) 
(lib) 



At p = this coincides with the previous results. Differ- 
entiating with respect to p and then setting p = Q gives 
a new result: 



dEr, 



dp 

dip 

dp 



dEn 



p=0 



p=0 



dp 
dp 



p=0 



p=0 



(12a) 
(12b) 



(There is a caveat that goes with these derivative rela- 
tions: It is possible that the derivatives do not exist. For 
example, one cannot rule out the existence of a branch 
point at p = 0, though there is no reas on to expect it.) 
Peshier and Thoma have obtained Eq. ( |l2a|) nonpertur- 
batively using a Ward identity argument ]40|] . In pertur- 
bative calculations the opposite slopes of the real energies 
is displayed in the one-loop calculations [l|J^-p|,|39t and 
the opposite slopes of the damping rates is a property of 
the calculations of Pisarski M and of Boyanovsky et al 



p3[ if one extends their answers to p = 0. Strictly speak- 
ing, the hard thermal loop calculations of the damping 
rates do not allow zero momentum but are limited to 
p > g^T. Repeated differentiation of Eq. (|ll] ) gives the 
general relations 



dp" 



p=0 



p=0 



y > Qpa 

^ ' dp" 



p=0 



p=0 



(13a) 
(13b) 



provided the derivatives exist. 

The vanishing of at p^ — £p and of Z3_ at po — 
along with the reflection symmetry (|^) means that the 
retarded propagator has the form 



S^'iPo^p) 



-J(7o-7-p)(^^-d(po,p) 
2 \pa-tp 

+ l{lo+l-p)(^^-d{po,^p)). (14) 

2 \pn-th ' 



Here d is an unknown function except for the requirement 
that it has no singularities in the upper half of the com- 
plex po plane. The complex residues satisfy the reflection 
property 



Zp{-p) = Zh{p), 



which implies 



Zp(0) = Z^(0) 
dp"- p=o dp" 



p=0 



(15) 

(16a) 
(16b) 



Equality of the residues and of the first derivatives has 
been found in one- loop perturbation theory 



6. Charge conjugation invariance 

The propagator in Eq. ( [1^ ) is invariant under chirality 
and parity by construction. Appendix C shows that in- 
variance under time reversal requires Eq. ( |C12| ), which is 
automatically satisfied by the ab ove. However invariance 
under charge conjugation, Eq. ( C13| ), is not automatic 
but requires the denominators to satisfy 



D.ipo,p) = -[D+{-p*,p)] 
Imposing this on Eq. (^4[) gives 

s^'ipo,?) 

1/ - -n/ Zp{p) 

= :^(7o-7-p)( c 
2 \po - tp 

+ ^{lQ+7-p)[ 



(17) 



PO + £*H 

z;{p) 



Po 



- fiPa,p)^ 

ri-Po,p)]. (18) 



Po - £h 
at i 

tiparticles of the hole excitation and particle excitation. 



The new po poles at —S^ and at —£* are due to the an- 
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respectively. Because £p and £h are in the fourth quad- 
rant, then —£* and — £^ are in the third quadrant. The 
function f{po,p) can have no singularities in the upper- 
half of the complex plane. It must satisfy the reflection 
property 



/(Po,p) = ~f*{~Pa,-p), 



(19) 



which is observed in the one-loop calculations 

The matrix structure of the propagator comes from the 
usual Dirac spinors: 



1/ 

o (70 - 7 • 



+ 7 ■p)^f^ = ^Va (-P, s)vp {-p,s). (20b) 



Both sets of spinors are eigenstates of x = 75^-^ = "foj-p. 
Since those in Eq. ( ^Oa| ) have x = +1, the chirality and 
helicity have the same si gn fo r particles and for anti- 
holes. The spinors in Eq. ( ^Oh| ) have x = — 1, indicating 
that the chirality and helicity have the opposite sign for 
antiparticles and for holes. 



7. Additional properties 

The adva nced propagator is obtained by the relation 
in Eq. (C5) of Appendix C: 



S^{po,P) 



= ^(70-7 

+ i(7o+7-p)(^ 
2 \po-t, 



z;{p) , z„{p) 



Pa + £h 
Zp{p) 



Po 



Sri 



f{~Po,p)]. (21) 



It has four poles in the upper half-plane as well as branch 
cuts in the upper half-plane from /. 

As discussed in Appendix C, the spectral function may 
be expressed as the difference between the retarded and 
advanced propagators. It is convenient to define that 
part of the spectral function proportional to 70 as 

PoiPo,p) = ^Tr[p(po,p)7o]- 
As shown in Appendix C, po must be positive: 
-Zp{p)2£p Zh{p)2£h' 



PoiPo,p) = -Im 
+Im 



Po 



f2 



Po 

f{Pa,p) + f{-Po,p) 



>0, 



(22) 



The canonical anticommutation relations of the fermion 
field operators impose require the integral of po to be 
unity as shown in Eq. (|C9|). This gives 



Re{Zp{p)) +Re{Zhip)) 



dpQ 



It is perhaps worth noting that the explicit one-loop cal- 
culations satisfy additional sum rules for the first and 
second moments of the energy, but these apply only at 
one- loop 

The time-ordered propagator can be expressed directly 
in terms of the retarded and advanced propagators as 

5"(po,p) = S"ipo,p)e''P«nipo) + S^{po,p)nipo). (24) 

It is rather trivial to prove this for po real. Appendix 
D proves that it holds throughout the complex po plane. 
Therefore in addition to the kinematic poles coming from 
the Fermi-Dirac function, the time-ordered propagator 
has eight dynamical poles: four in the lower half-plane 
from and four in the upper half-plane from S^. Of 
course the dynamical poles are all reflections of the basic 
particle and hole energies £p and £h- 



8. Conclusion 

The above results depend only upon invariance un- 
der chirality, parity, charge conjugation, and time rever- 
sal and not at all upon perturbation theory. The com- 
plex energies £p{p) and Ship) the particle and hole 
excitations are gauge-fixing invariant as proven generally 
by Kobes, Kunsattter, and Rebhan However the 

residue functions Zp and Z^ are expected to change with 
gauge. Appendix D shows that the renormalized electric 
charge for the particle and hole excitations, as measured 
by the coupling to photons at zero momentum, has the 
correct value independent of the functions Zp and Zh- 

The general shape of the two dispersion curves requires 
difficult calculation. The only simple property is that 
as p — > 00 the effects of temperature diminish so that 
£p{p) — !■ p and Zh{p) — > 0. The hard thermal loop calcu- 
lations enjoy several properties that have not been proven 
to hold generally. Namely, in the HTL approximation 
the hole energy is asymptotic to p as p — > 00, the phase 
velocities of both excitations are larger than unity, and 
the group velocities of both excitations are smaller than 
unity. 

At temperatures close to the critical temperature the 
running coupling is large and one-loop calculations do 
not apply. Since the doubled dispersion relation is such 
a characteristic feature of the high-temperature, chirally 
symmetric phase of QCD, it is frustrating that experi- 
mental signatures are so difficult to find. One possibil- 
ity are the Van Hove singularities in dilepton production 
found by Braaten, Pisarski, and Yuan [O. This has 



Im/(po,p) = l. (23) 



been recently studied by Schafer and Thoma |39| and 
Peshier and Thoma |^^. The Van Hove singularities are 
determined by the minimum value of the hole dispersion 
relation i?„iin and by the maximum energy difference be- 
tween the two dispersion relations AE. In the dilep- 
ton rest frame the dilepton rate has a square root diver- 
gence at /co — 2£'min and at — Ai?. Unfortunately in 
the HTL approximation the values -Bmin — 0.38gT and 
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/S.E = 0.19gT are so small that the effect is swamped by 
the large dilepton continuum. However since the contin- 
uum falls rapidly with energy, it is possible that the Van 
Hove singularities might exceed the continuum if the true 
values of i?min and E'max are large enough. Of course, the 
effect will be weakened by the quark damping rates. 
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APPENDIX A: p = independent of p 

A crucial step in the proof that particles and holes have 
separate dispersion relations, which is used in Eq. and 
subsequently, is the fact that at zero thrcc-momcntum 
the fermion propagator does not depend on the momen- 
tum direction p. (This property is specific to fermions. 
Gauge boson propagators can depend on unit vectors at 
zero momentum, as occurs even in the free Coulomb- 
gauge propagator.) The argument is simple. The self- 
energy at any order of perturbation theory is expressed 
in terms of momentum integrations over integrands com- 
posed of fermion propagators and boson propagators. 
The routing of the external three momenta p through 
the internal propagators depends on what choice is made 
for the loop momenta kj that must be integrated. Be- 
cause of fermion number conservation at each vertex, the 
external fermion number can be uniquely traced through 
a sequence of internal fermion propagators that form a 
continuous path through any self-energy diagram. There- 
fore one can always choose the loop momenta so that the 
gauge boson propagators and the fermion propagators in 
closed loops will contain various kj but never p. The ar- 
gument of the linked fermion propagators that connect 
with the external lines will be the sum of p with a linear 
combination of the kj . These internal fermion propaga- 
tors depend on the vector p and not separately on p and 
\p\. Hence if the external momenta \p\ is set equal to 
zero, there can be no dependence on the direction p and 
the propagator must be as shown in Eq. (p|). 



APPENDIX B: Symmetry under p 

It is convenient to parameterize p by p, 9, c 

Px = P sin 9 cos (j) 
Py = p sin 9 sin (j) 
Pz = pcos9. 



Under the transformation p— > —p, 0— s-vr — 0, 0— !-(/>-f7r 
the vector p is unchanged. This amounts to p ~p and 
p —p. If the propagator is considered as a function of 
the four variables po,p, and p, it must not change: 

S"{po,P,p) = S^{po,~p,~p). 
In terms of the definitions in Eq. this requires 

D^ipo,p)^D+ipo,-p), (Bl) 
as employed in Eq. (^). 

APPENDIX C: Retarded propagator 

A natural starting point for discussion of any finite- 
temperature propagator is the finite-temperature spec- 
tral function: 



Pa/six) =^e~ 



I3E„ 



{n\{Mx),i,p{Q)}\n) 
Tr[e"/3«] 



(CI) 



where the sum is over a complete set of energy eigenstates 
\n). In the following, Dirac indices will be suppressed and 
p will be treated as a matrix. The matrix satisfies: 

[p(a;)]t=7op(-x)7o. (C2) 



1. Retarded vs advanced 

The basic retarded and advanced propagators are de- 
fined in coordinate space as 

S^{x) = -i9{xQ)p{x) 
S'^ix) = i9i-xo)p{x). 

In momentum space they have the simple dispersion re- 
lations 



s-iPo,Pi=r 

J-oo 

S^{P^,P)-- 



duj p[uj,p) 



po- uj + ir] 
duj p([x),p) 
2ti po ^ uj — ir] 



(C3) 
(C4) 



The retarded propagator is an analytic function of com- 
plex pq in the half-plane Impo > 0. Its only singulari- 
ties are in the lower half-plane. Similarly the advanced 
propagator is an analytic function of complex pq in the 
half-plane Impo < 0. All its singularities are in the up- 
per half-plane. The two propagators are connected by 
the relation 



5^(pS,p)=7o[5«(Po,p)]^7o, 



(C5) 



which follows directly fr om the adjoint property of the 
spectral function in Eq. (C2). 
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2. Spectral Function 



Using Eqns. (|C3|) and ( p4D the spectral function in 
momentum space can be expressed as 



p{iu,p)^t[S''iiu,p)-S^{uo,p)\. 



(C6) 



Although the spectral function is defined for real uj, this 
relation allows uj to be continued into the complex plane. 
From the definition of the spectral function the ma- 
trix [/5(w,p)7o]q;/3 is Hermitian and has a positive-definite 
trace: 



po{uj,p) = ^Tt[p{uj,p)jo] > 0. 



(C7) 



This leads to the condition in Eq. (^2|). 

From the canonical anticommutation relations at equal 
times, the spectral function must satisfy the relation 
jqS^{x) = p(0,x). The Fourier transform is 



7o 



duj 
2^ 



p{uj,p}. 



From Eq. (C6) this can be expressed as 



70 = « 



(C8) 



When Eq. and (|2l| ) are used the terms proportional 
to 7 • p automatically integrate to zero and those propor- 
tional to 7o give the sum rule 

/OO 7 
^Po{^,p)^l. (C9) 



3. Discrete symmetries 

In the high temperature phase, QCD is expected to be 
invariant under chirality as well as under parity, charge 
conjugation, and time reversal. These symmetries con- 
strain the spectral function and consequently the propa- 
gators. 

(i) Chirality invariance requires that the spectral func- 
tion satisfy p{x) = —J5p{x)^5. For the retarded propa- 
gator this implies that 

Qs: S''{po,p)^-j5S''{po,P)j5- (CIO) 

(ii) Parity requires p{t, x) ~ jop{t, ~x)jo and therefore 
P : S^{po,p) = 7o5^(po, -P)10- (Cll) 

(iii) Time reversal invariance requires [p{t,x)]* = 
Tp{—t,x)T where T = 17-^7^ and relates the retarded 
propagator to the advanced: 

[S^iPo,P)r =TS^ip*o,-p)T. 



Using Eq. ( |C5| ) this can be expressed in terms of S^: 
T : 5«(po,p) = T7o[5«(po, ~P)floT. (C12) 

This relation is automatically satisfied by Eq. ([T^). 

(iv) Charge conjugation requires = Cp{—x)C 

where C = 17^7'^. This relates the retarded propagator 
to the advanced: 

[S''{po,P)f = -CS^{~po,~p)C. 



However, using Eq. ( [Cq ) this can be rewritten as a con- 
straint on the retarded propagator directly: 

C : S^{po,p) = -C7o[5^(-pS, -p)]*loC. (C13) 

It is this relation that leads from Eq. ( jl^ ) to Eq. 

(v) It is sometimes convenient to deal with the com- 
bination 6 = TCP, which is always an invariance. Thus 
[^(2:)]^^ — ~7o75P(2;)757o always holds and therefore 

TCP: [5«(po,p)]^= 757o^''(-pS,-p)7o75 (C14) 
must always hold. 

APPENDIX D: Time-ordered propagator 

This appendix proves Eq. (|2^), which expresses the 
time-ordered propagator in momentum space directly in 
terms of the retarded and advanced propagators. The 
relation holds for all po . The time-ordered propagator is 
defined as 



(n|r(^„(a;),V0(O))|n) 



Tr[e-'3^] 

This can be expressed in terms of the spectral function 
cLu / p{uj , p) [I ~ n{LL!)] p{uj,p)n{u)) 



S''{po,P} 



-00 27r V P0 + iij - oj 



Pa - i-q-u) 



where n(w) is the Fermi-Dirac distribution function: 



n{Lj) 



1 



exp(/3w) 



It is rather easy to evaluate S^^ for pg real and 77 in- 
finitesimal and obtain Eq. (|2j) for real po- However, 
since the focus of the present paper is the pole structure 
of the propagators in the complex pg plane, that easy 
argument is not adequate. In order for Eq. (U) to hold 
in the complex po plane it will be important that 
have not branch points on the real axis. 

The dispersion relation for S^^ shows it to be the sum 
of two functions. The first is analytic in Impo > and 
the second is analytic in Impo < 0. Write S^^ as 

^"(po,p) - 5^(po,p) + F(po,p) + G(po,p), 
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where 

Impo > : F{po,p) 



— oo 
oo 



dio p{uj,p)n{uj) 
27r po — Lo 
duj p{uj,p) n{uj) 



2-K 



Po 



Using Eqns. (C6) for the spectral function, F can be 
computed by closing the uj contour in the upper half- 
plane for the term containing and closing it in the 
lower half-plane for the term containing S^. Since n^uj) 
has poles at uj = ztfte where ili = iiirT the result for F 



^ t , rr.^ f S''{ne,p) , s^{-ni,p} 

-S'^{pa,p)n{pa) - T^[ — + — — 



Po 



This is analytic for Impo > but has poles and cuts for 
Impo < 0. In the same fashion G can be computed by 
closing the lo contour in the upper half-plane for the term 
containing and closing in the lower half-plane for the 
term containing 5'"^ with the result 



T 



s^{ni,p) , s^{-n,,p) 



V nn- 



Po -^e Po + 



+ S^{pa,p}n{pQ). 



This is analytic for Impo < but has poles and cuts in 
Impo > 0. When the results for F and G are added the 
sums over £ coming from the poles of the Fermi-Dirac 
function n{uj) cancel and give 

F{po,p) + G{po,p} = ( - 5^(po,p) + ^^(Po,R)) n{po). 

The final result for the time-ordered propagator is 

S^\po,p) = 5^(po,p) [1 - n{po)] + S^{po,p) n{po). 

Since 1 — n{po) = e^'P°n{po) the result can be expressed 
asEq. (H). 

APPENDIX E: Electric charge 

This appendix will show that the four poles in the re- 
tarded propagator (jl^) all have the correct value of the 
renormalized electric charge. The operator ^ for the u 
quark destroys charge = 2e/3; the operator ij) for the 
d quark destroys charge = — e/3. The amplitude for a 
photon with zero four-momentum to couple to a fermion 
is 

Sr{pq,p)T ^,(jpq,p)Sr{pq,p). 
At finite temperature the Ward identity 

determines the photon vertex function Q-Q. The in- 
verse of the propagator in Eq. (S) can be written 



[S{po,P)\ ^ = loA{Po,p) - 7 ■ pB{po,p), 
where D± ~ B. The Tq vertex is thus 

ro(po,p) ^ joA'{po,p) - J ■ pB'{pQ,p), (El) 
where prime denotes the po derivative. 

1. Pole at £p or 

When Po £p the retarded propagator, Eq. (|l8|), is 
dominated by the pole 



S^{po,p) Zp{p) 



,u{p)u{p) 



Po - Sp 



If the quark field "0 carries bare electric charge eo^ , then 
the renormalized electric charge which couples the pho- 
ton to the particle excitation is 

Q = eQ^^/z^Zp{p) ■ u{p}ro{pQ,p)u{p), 

where \^Z^ comes from the photon propagator and ^/Z^ 
comes from each of the two fermion fields at the vertex. 



Substituting Eq. (|El|) gives 

Q = eo^y^Zp{p)[A'{pQ,p) - B'{po,p)]^^^^^. 

As Pa — > £p the denominator = A ^ B vanishes lin- 
early: A — B ^ {po — £p)IZp[p)i so that 



It is easy to check that the pole in Eq. 
also couples with the same charge (E2) 



at Po 



(E2) 



2. Pole at £h or —£p 



At the hole excitation, pq £h, the retarded propa- 
gator, Eq. (|l8|), behaves as 



S"ipo,p} ^ Zh{p) 



J2spin'"i-P)'"i-P) 



Po - 

The renormalized electric charge of the hole excitation is 

Q = eQ-^y/zlZh{p) ■ v{-p}Tn{po,p)v{-p). 
Substituting from Eq. ( pi] ) gives 

Q = eQ^\f%,Zh{p) \A!{po,p) -HB'(po,p)]pj,^£:,^- 

At the hole excitation, L)_ = A + B ^ {pq — Eh)/Zh{p) 
as Pq — > £h so that 



(E3) 



The same argument applies to the pole in Eq. (jT^) at 



Po 



-£■1 
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